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A  Computer  Subroutine  for  Stress 
-'J' 
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by  John  E.  Brock 

Based  upon  theory  developed  by  the  -writer,  R.  E.  Brown  developed 
a  successful  computer  program  for  analysis  of  radial  and  circumfer- 
ential stresses  in  rotating  axisymmetric  disks  of  variable  thickness 
having  an  axisymmetricai  thermal  strain  field.  The  writer  revised 
Brown's  program  so  as  to  invoke  a  group  of  ancillary  subroutines 
which  have  been  found  useful  in  another  application.  In  doing  so, 
however,  much  unnecessary  and  confusing  normalization  was  introduced. 
In  particular,  one  of  the  normalizations  would  cause  the  analysis  to 
fail  in  the  quite  common  case  of  a  disk  with  no  radial  loading  at  its 
outer  boundary.  All  this  material  appears  as  Reference  i,  hereof. 

Referees  evaluating  a  paper  based  upon  Reference  I,  called 
attention  to  these  faults  so  that  the  program  has  been  rewritten. 
A  listing  of  the  main  subroutine,  RODISX,  as  revised,  as  well  as 
listings  of  the  ancillary  subroutines  may  be  found  in  Appendix  A 
hereof.  The  reader  will  note  that  other  changes  have  also  been  made 
resulting  in  somewhat  more  flexibility  of  application.  Employment  of 
the  revised  program  is  described  in  the  textual  material  which  appears 
at  the  beginning  of  the  listing. 

Appendix  3  contains  a  revision  of  the  second  illustrative  example 
problem  of  Reference  I.  This  problem  was  solved  for  various  values  of 
M  =  N-l,  the  number  of  equal  subdivisions  into  which  the  annular  radius 
b-a  is  divided  for  purposes  of  numerical  analysis  by  R0DI3K.  Also,  a 


number  of  different  values  of  KP(3)  were  used.  If  K?(3)  >  0,  its 
value  is  the  number  of  iterations  which  will  be  performed  by  RODISK. 
If  XP(3)  <  0,  iteration  will  continue  until  three  successive  values 
of  che  unknown  parameter  3  determined  in  the  course  of  the  analysis, 
satisfy  the  relation 

|Bi-32|  +  |32-33 !  +  [B3-Bi|       ,^m 


Bi   +  B2  +  B3 


We  also  determined  execution  time  by  use  of  the  library  subroutine 

IXCLOK,  executing  under  CP-cms  on  the  IBM  360/67  at  the  W,  R.  Church 

Computer  Center  at  the  Naval  Postgraduate  School. 

We  found  that  execution  time  per  iteration  is 

t_   =  1.2  M  +  5   (milliseconds) 
iter 

for  any  problem. 

Accuracy  was  evaluated  by  dealing  with  problems  raving  available 

analytic  solutions.  It  was  found  that  the  principal  limitation  on 

accuracy  is  determined  by  the  choice  of  subdivions,  the  integer  M   = 

M-l,  so  that  there  is  a  certain  inherent  error  regardless  of  how 

many  iterations  are  made.  This  error  depends  on  M,  of  course,  and  upon 

the  details  of  the  problem.  The  error  is  greatest  near  the  inner 

radius  of  an  annular  disk,  and  is  large  if  the  ratio  a/b  is  small. 

Fortuitously,  the  error  may  be  smaller  for  an  early  iteration  than 

for  a  somewhat  later  iteration  but  this  is  not  practically  useful 

information.  For  the  problem  of  Appendix  3  hereof,  with  a/b  =  .165, 

we  find  the  results  given  in  Table  1,  (see  next  page). 

Thus,  for  example,  with  M  =  20,  there  is  an  inherent  error  of 
about  1%   and  the  results  are  not  significantly  improved  by  iterating 
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approx . 

approx. 

total 

n 

limiting 
%   error 

iters. 
ren '  ci 

time, 
sees. 

5 

15 

5 

.055 

10 

5 

7 

.12 

20 

1 

11 

.32 

40 

.1 

17 

.90 

100 

.01 

25 

3.1 

Table  1.  Percent  error,  required  iterations,  and 
execution  time  for  problem  of  Appendix  3. 

more  than  eleven  times.  With  eleven  iterations,  the  solution 
is  returned  from  RODISK  in  0.32  seconds. 

The  significant  conclusion  is  that  the  execution  is  so  fast 
that  one  may  as  well  take  M  =  100  (corresponding  to  N  =  101,  the 
maximum  available  under  present  dimensioning)  and  iterate  many 
more  times  than  is  strictly  necessary.   Taking  N  =  101  and 
KP(3)  =  -3  gave  execution  in  3-7  seconds  with  31  iterations  and 
with  an  accuracy  of  0.00^$    (In  the  problem  at  hand,  a   (a)  was 
specified  as  zero  and  the  program  gets  -1.14E-11  so  that  the  error 
here  is  "infinite'.'  Our  evaluation  of  0.00^%   is  for  the  first 
position  rather  than  for  the  zeroeth.) 

This  concludes  the  text  proper  of  the  present  report .  However, 
we  take  advantage  of  this  opportunity  to  correct  errors  in  Reference 
1 ,  viz . : 

(1)  Page  3,  equation  12  should  read 

m  =  ±/(n2-4vn+4)  =  ±/[(n-2)2+4(l-v)n] 
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(2)  Page  6,  line  2.  In  place  of  T  read  oT. 

(3)  Page  6,  equation  33.  Lower  limit  of  integration  should 
be  a  rather  than  0 . 

(4)  Page  7,  line  following  equation  40,  Reference  should  be 
to  equation  37  rather  than  equation  33. 

Acknowledgment  is  gratefully  made  for  assistance  by  the  Naval  Post- 
graduate School  Research  Foundation.  Appreciation  is  also  expressed 
to  the  referees  of  the  ASME  Journal  of  Applied  Mechanics  for  directing 
attention  to  the  flaws  in  the  earlier  version  of  RODISK. 
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Appendix  A 

Listing  of 
Subroutine 
RODISK 
and  ancillary 
subroutines 


C  SUBROUTINE    RODISK.       JOHN     E.     EROCK,     1     MAY    1978,     REVISED    1    AUGUST    1<57£.     FCDCCC1C 

C  THIS    IS    A    SUePCUTINE    r  CR    DETERMINING    RJOIAl.    ANO    C  IRCU  NFEREN  TI  AL     STR ESS  ROC  COG  20 

C  IN    AN    AXI  SYMMETRIC    THIN     ELASTIC    DISK     HAVING     AN     AXISYMMSTRIC    THERMAL  ROOCCC2C 

C  STRAIN    FIELD   AND    ROTATING    AT    ANGULAR    VELOCITY     ONEGA     (R4DIANS/SEC0NDJ        RCDCCCAO 

C  ABOUT    THE     AXIS    OF    SYMMETRY.       TWO    TYPES     OF     PROBLEM    MAY     BE    TREATED:  RCC0005G 

C  TYPE    U     ANNULAR    DISK    CF     IN  SI  CE    R A C I  US    ARAC    AND     OUTSIDE     RADIUS  RCCCCCtC 

C  BRAD.        THE    RACIAL     STRESS     IS     ERA    AT    THE     INNER    RADIUS    AND  RCCCCC7C 

C  SR3    AT    THE    OUTER    RADIUS.       THE     INSICE    RACIUS    MUST    BE  RCOCCCcC 

C  GREATER    THAN    ZERO  R0CG0C9G 

C  T#FE    2:     SCLIC    OISK    HAVING   RADIAL    STRESS    SRB    AT  OUTSIDE    RADIUS    BRAD.     RCOOOIOD 

C  THE    USER    MUST    PROVIDE     A     "AIN    PROGRAM     VHICh    CALLS     SUBROUTINE    RODISK  RCDOC110 

C  AFTER     IT    HAS    SUPPLIED    THE     fCLLCWING     IN  FORMAT     ION.  PCCCC12G 

C  (1)     Nf     INTEGER.     (  N- 1  )     IS     THE    NUMBER    OF    EQUAL    SUBDIVISIONS    INTO    Wh  I CHP  CC  CO  1 20 

C  THE     ANNULAR    RADIUS     (3PAO    MINUS     ARAL")      IS     OIVICEC     FCP    CCMPU-  ROCCGl^O 

C  TATIONAL    FURPCSES.     THE    PRESENT    DIMENSIONING    CAtf    ACCOMMODATE    N  RCCCC1SC 

C  NOT    GREATER     THAN    101.  ROD00160 

C  (2)     BRAD  RCDCC170 

C  (3)     ARAC     (NOT    NECESSARY     FOR    PROBLEMS    OF     TYPE     2.)  RCCCC16C 

C  U>     SRB  RODCC1S0 

C  (5)     SRA     (NOT    NECESSARY     FOR    PROBLEMS    QF    TYPE     2.)  RCDCC2C0 

C  (6)     POIS,     FOISSON'S    RATIO  RODC0213 

C  (J)     KP(1)=1,2.       INTEGER    TO    DEMOTE    PROBLEM    OF    TYPE    1,2.  RGDCC22C 

C  (8)    KP(2),     INTEGER    TO    PROVIDE    FOR     SKIPPING     kHILE    PRINTING  RCCCC22C 

C  OUTPUT.       FCR    EXAMPLE,     IF    N=101     AND    KF(2)  =  5,     ONLY     EVERY  PCD  C  C  2^0 

C  FIFTH    SET    OF     VALUES-     kl  LL    BE     PRINTED:    1ST,6TH,...,     96TH,  RCC002fO 

C  AND    101ST.  FGDCC260 

C  (9.)     KP(3),     INTEGER    S  PEC  I F  YING    NUM  E  ER    OF     ITERATIONS    TOSE  RCCCC2TC 

C  PERFORMED.       USUALLY     *P!3  i  =  10     IS     SUFFICIENT     FGR    SNGIN-  ROCCC280 

C  EERING    ACCURACY.       ALTERNATELY,     IF     Kp(3)     IS    A     NEGATIVE  ROCCQ2<=0 

C  INTEGER,     ITERATION     WILL    CONTINUE    UNTIL    THREE    SUCCESSIVE  PCDCCECC 

C  VALUES    OF    A    PARAMETER     E,     OETERVINEC     INTERNALLY,    ARE  RCECCE1C 

C  SUFFICIENTLY    CLCSE     ^     COMPARED     TO    AN    ERSlLCN     EQUAL     TO                      .  FCDCC22C 

C  TEN    RAISED    TO    THE     KP<2)     POwER.  RCCCC22C 

C  (1X3)     KP(4).     IF    ,<P(4)=0,     ONLY     FINAL    ANSWERS     WILL     EE     PRINTED.  RCDCC24C 

C  IF     KP(4)=1,     A    SEQUENCE    OF    ITERANT     VALUES     QF     E    WILL     BE  RCCCC25C 

C  PRINTED    TO    INCICATE    DEGREE    CF    CONVERGENCE.       IF    K  P  (  4  )  >  1  ,  FCCCC36C 

C  THERE     WILL    BE    NO    PRINTING    AT     ALL    WITHIN    RODISK,     SUT    UPON  PCCCC27C 

C  RETURN    KF(5)     WILL     CONTAIN    THE     NUNEER     OF     ITERATIONS     VHICH  RCCCC2EC 

C  WERE     FERFGRMEO    SO    THAT    KP  (  5  )     MUST    BE    RESET    BEFORE    RODISK  RCDCC2?C 

C  IS    CALLEO    AGAIN.  RCDCQ4CQ 

C  (11)     KP(5).       KP(5»=0     CALEBS    MRNTS     CUBIC     SPLINE     INTEGRATION  RCCCC410 

C  TO    BE    USED.    OTHERWISE    TRAPEZOICAL     INTEGRATION     IS    USED.  RCDCC42C 

C  (1.2J     VECTOR    X(1,J),     J=1,2,-..,N,     CONTAINING     VALUES    OF     DISK  FCCCC42C 

C  THICKNESS    AT    EQUALLY     SPACED    RACII     FR.CN     INSICE     TO    OUTSIDE.  RCDCC44C 

C  (13)     VECTOR     X(2,J)     CONTAINS    VALUES    CF    GAMMA     TIMES    ONEGA  FCCCG4SC 

C  SQUARED    WHERE     GAMMA     IS     (MASS)     DENSITY     Oc     THE    MATERIAL.  PCOCCAdC 

C  FOR    MOST     PROBLEMS     GAMMA    OGES    NOT    VARY     WITH    RADIUS    AND  FCCCC47C 

C  ALL    ELEMENTS    OF     THE     VECTOR     WILL     BE     THE     SANE.  PCDCCAtC 

C  (14)     VECTOR     X(2,J)     CONTAINS    VALUES     CF     (  EE  )  ( A L F HA )  (T E E  )     WHERE  PCCCC4=C 

C  EE     IS     YOUNG'S    MODULUS,     ALPHA     IS     THE    COEFFICIENT     OF     LINEAR  RCCCC5CC 

C  THERMAL     EXPANSION,     it^D    t«     is     THE     TEMPERATURE    CHANGE.  RCOCCS1C 

C  THE    MAIN    PROGRAM.    MUST    CONTAIN    THE    STATEMENTS:  PCOCC520 

C  IMPLICIT    PEAL*8     (A-H,C-Z)                                                                                                      •     PCDCC33C 
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WRIT; (6,  2)  KPil) 
RQOISK  PROBLEM  OF 


FEAL*8    X  (  20,  101 ) 

INTEGER    KP(5) 

COMMON    X,N,KP 

CCMMCN    /CNE/ARAQ.ERAO,SRA,SR3,P0!< 

FOLLCWING    SUBROUTINE    RGCI3K    THERE    ARE     SEVERAL    ANCILLARY 
SUBROUTINES    WHICH    PERFORM    vARICUS    CPSFATICNS    CN    THc    VICTORS 
X(1,J).        THE     PURPOSE    OF     E/CH     13     GBVICUS    FRCN     THE     LISTING. 
THEY  MAY    BE     EMPLOYED    IN     THE    USER'S    Mil  I N    PRCCRAM.        SUEPCUTINE 
OUPV,     WHICH    CUPL1CATES    A     VECTOR,     ANO     SU8RCUUNE     PRIV,     «HICH 
PRINTS    A     VECTOR,    ARE     WC7    CALLED     EY     PODISK     BLT    MAY    6c     USEFUL 
IN     THE    USER'S    MAIN    PROGRAM. 

SL8RQUTINE    RC3ISK 

I  fFLICI  T    REAL*d     (A-HtQ-2  J 

REAL*8     X(20,1C1) 

INTEGER    KP(5) 

CCMMCN     X,N,KP 

CCMMGN     /CN  E/ ARACSRAQ,  SPA, SR3, POIS 

0N£  =  1  .D-»0 

ZERO=O.C-fO 

81=1  .0  +  3 

B2-1.0  +  10 

33=1.0+12 

1F(KP(3).LT.0)     EFS=(1.C+I)**(KP(3)) 

I  F(  KP<  1)  .EQ.  2)     ARAG=ZERC 

EMA=BRAO-ARAO 

ENM=N-1 

I  f  UP  (<,]  ,LE.l  ) 
2    FCRMATt//  ,  10X  ,'ROuISK     PROBLEM    OF     TYPE     ',11,'.  «) 

CO    5     1=1 ,N 

E1M=  I-  1 

Y=EIM/ENM 

X  (4,1  J  =AR AC + (8R AD -A fi AC)* Y 

X(5tl J =Y 
5    X  (6,1  )  =Y 

I  TE  R=  1 

I  F<  KP(  1  )  .£0.2)    GO    TO     ICC 
THE     FFGBLEM     IS     GF   TYPE    1:     ANNULAR    DISK 

Cl=(2.D+0+P0IS)*<8RAC-ARAD) 

CALL     INTV(  1,7  ,8MA) 

C2=X(7,N) 

C5=X(2,  ])-X(3  ,N)  +  (CNE-FOIS)*(SRA-SRE) 

CALL    MUL V<1 ,2,8) 

CALL    MULV( 4,8,9 ) 

CALL     INTV(9,10,2MA) 

C£  =  X(  10  ,N>+X(i,N)*SRB-XU,l)*SRA 
20    CALL    I NTV ( i  ,11, BMA) 

O8RA0+ (CM E+PQ  I  S)  *X(  1  1,N  ) 

CALL    MU  LV(   1,  6,  12) 

CALL    I NTV(12,  13, SMA) 

C4  =  X( 12, N) 

D=C1*C4-C2*C3 

A=<C5*C4-C6*C3) /O 

B*(C1*C6-C2*C5)  /D 
20    CCNTINUE 

IF  (KP(A).EC.l)     WRITER, T)     ITER,  A,  8 
7     FCRMAT( 5>, I10,1F2£20. 5) 

CALL    MUCS(7,  14, A) 

CALL    MULS(13, 15, E) 

CALL    AOCV(  14,  15,  15) 

CALL    SUBV( 15,  10,  16) 

S=SR8*X(  1, N)-X( 16 ,N  J 

CALL    AODSi  U,lc,S) 

CALL    OIVV(  1c,  1,  16) 

ZA  =  X  (3  ,  l)  +  A*ARAQ+(CNE'-fJ0IS  )  *X  (16,  1) 

CALL    MULS(11,17,8 ) 

CALL    SUES ( A  ,13 ,ARAC ) 

CALL    MULSf 13,13 ,A ) 

CALL    AOOV(  17,13,1  1) 

$  =  -< CNE+FOIS  ) 

CALL    MULS(  17,  17, S  ) 

CALL    ADDS  (17, 17 ,ZA) 

S=CNE-PCIS 

CALL    MULS (16 ,13, Si 

CALL    SU8V(  17,3,  17) 


FCCCC540 
RCOCC55C 
RCC00560 
FO0CC5  70 
RCCCCS6C 
FGDCCS^C 
RCDCC6CC 
ROC006  10 
F CDC 06  20 

RocccdEc 

ROCC0640 
ROCCC65C 
RO  CCC660 
RCCCC670 
R0CCC6E0 
RCC00690 
R0CCC7CC 
RCDCC71C 
R0CC0720 
R CD  COT  30 
RGC  C  C  "i  AC 
ROCCC750 
RCCC076C 
RO  CCC77C 
RCDC0780 
R00CC7SC 
RGOCCcCO 
ROCCCE  1C 
RCDCCE2C 
RCCC082J 
FCDG0840 
RGCCC65C 
P0CCCE6C 
RCCCCtTC 
ROCCC65C 
RCCCCE^C 
ROC  CCS CO 
RCCCC910 
RCCCC92C 
RCDCC=2C 
RGDC09-<0 
FCDC0950 
RCCCC^C 
P00CC970 
RGDCC«iC 
ROCCC^O 
PCDC10C0 
ROC  CI  CI C 
R0CC1C2Q 
RGCC1C2C 
R0DC1C4C 
RCCC1C5G 
PODC1C60 
RCDC1CTC 
P0CC1CEQ 
RCDCIC'C 
R0CC11CQ 
PCDC1110 
RCCC112C 
PCDC112C 
R0CC114C 
RO CC11  50 
PCCC1160 
RO001170 
RCDC11EQ 
ROCC  11=0 
RCDC12C0 
RCCC1210 
FCOC1220 
RC0C122C 
PCOC  12^0 
RCCC125C 
ROCC1260 
P0CC1270 


-     6     - 


C/LL    SUEV(  17  ,  18  ,17 )  P0CC12S0 

I7ER=ITER+1  RCDC12<C 

I  F<KP(  3  )  .LT.O  )    GJ    TG    215  RQCC12C3 

IFU7ER.GT.XP  (3)1     GJ    TG     200  PCDC13  20 

40    CONTINUE  RO0C152C 

CALL     DUPV<4,14)  ROCC1330 

I  F(KP(  1)  .£0.2)    X(14,1)  =  C!\E  R0CC134C 

CALL    DI  VV  (  17,1  4,1  6)  ROCC1350 

CALL    SUBS  I ic,  IS, A)  RCCC1360 

CALL    DI  \.S  (13,  6,81  RODC1270 

1F(KP(  1  )  .EQ.2)    GO    TO     15C  RCDC13S0 

GC    TO    20  ROCC1290 

100    A=ZSRC  RCDC14CC 

C=(ONE-POIS  1/X<  1 , 1 J  R0D01410 

CALL    MULV(l,2i7)  PCD01420 

CALL    MULV(4,7,8)  RC0C142C 

CALL     INTV (3  ,  10 , 2MA )  FCCC1440 

SUM=X(10,N)+SR6*(X(1,N)-X(1,1))  RC0C14  5C 

Siy=C*  SUM+X  (2,1)-X(3,N,J  RO  CC146G 

150    CALL    I NTV( 6,11,3MA)  PCCC147C 

DEN  =  6RAO(0NE+P0IS)*X<il,N)  R0D014K 

CALL    MULV(1,6,11)  R0DQ1490 

C/LL     INTV(  11,  13,  3AA  )  RCCC15CC 

DEN=OEN+C*X  (13, N)  R0DC151G 

B  =  SUM/OE.N  PCCC1520 

CALL    !NTV(1,7,BMA>  F0DC1520 

CALL     INTV(6,11,BMA)  RCDC154C 

GC    TO    30  ROCClffG 

215    B3=B2  RCDC156C 

82=81  ROC01f70 

81=8  PC0C1530 

ZUM=OAeS(  e  1-82)  +0A8S  (  E2-E3) +  CA6S(  3  3-8  1  >  ROC  CI  5 CC 

DlV=OA3S(61)+0AeS(32)+CA3S(63)  P0D016C0 

CRIT=ZUM/CIV  R0CC161C 

I  FiCRI  T.LT.SPS)     GC    TO    200  R0DCU2G 

GO    TO    4  0  PCCC16  2C 

2CC    CALL    ACCV(17,16  ,19;  RODC1640 

I  f  UP(  3)  .LT.G.ANO  .KP(4  J.EQ.O)     WR  ITE ( 6 , 201 )     ITER, EPS  RGCC165Q 

201    FORMAT  l//t20X,I8,  •     ITERATIONS    RECUIREO    WITH    EPSILCN    =    ',1PE8.D  R0DC166G 

IF(KP(4  J.LE.l)     »RITS(6  t2'-4)  RCCCU7C 

204  FGRMATl//)  RCCCU60 
IF(KP(4)  .LE.l  )     WRITE(  6,  2C5)  FCDC1690 

205  FCRMAT( 23X,'RA0IUS'  ,10>,  'THICKNESS'  ,5X, • GAMNA    Cf-EGA     SQ«,7X,  RCCC17CG 
l'EE    ALPHA    TEE'  ,7X  ,«  SIGMA    RAD  IAL  •  ,  6X  ♦  •  SI  G.XA    CIRCUMF')  FODO  17  10 

NSKIP=KF(2)  RCCC172C 

00   210     I  =  1,N,NSKIP  ROCC173C 

J=I/NSKIF  RCCC174C 

IF<KP(4)  .LE.l)    WRITE(6,  211  )J  ,  XI  4t  1)  ,X(  1,1)  ,X<  2,1  )  ,  ROCC175G 

1X( 3,1  )  ,  X(16,l)  ,  X(19,I  )  PGCC176C 

210  CONTINUE  RG0C171C 
IF(KP(4)  .GT.l)     K?(5)=ITER  RCDC17EC 

211  FCRMATl  I10tlP6£19.5i  R0CC1790 
RETURN  •  PCDC18CG 
END  RCDC1S1C 

THIS    IS    THE    START    OF    THE    ANCILLARIES  ROCC1820 

SUBROUTINE    ADDV ( Nl ,N2 ,N3 J  RCDC163G 

REAL*8    X(20,10i),S  ROC016^0 

INTEGER    KP(5)  PCDC1850 

COMMON    X,N,KP  RC0C1J6C 

00    1    1  =  1, N  RO0C1870 

1    X  (N3,  I  )=X(N1,  IMX(N2,  I  )  RCCC1€£G 

RETURN  RODC1£90 

ENC  /  RCDC19C0 

SL6R0UTINE    SU8V  (  N  1  ,N2  ,N  :  )  PCDC1910 

REAL*8     ><20,101),S  RCDC1C.2C 

INTEGER    KP(5)  FCCC  192C 

COMMON    X,N,KP  RCDC1S4C 

00    1    1  =  1, N  R0CC19iG 

1    X  (N3,  I  )=X(Nl,i)-X  (N2,  I)  FCCC1960 

RETURN  RODCl^C 

ENO  RO0C19S0 

SUBROUTINE    MULV ( Nl ,N2 ,N2 )  P0CC19CC 

REAL* 3     X(2G,101),S  RODC2CC0 

INTEGER    KP(5)  RCD02010 
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CCMMGN    X,N,KP                                   •  ROCC2C20 

CC     I    1  =  It  N  S0CC2C3C 

X  <  N3  , 1  I  =  X  (  N  1 , 1  )  *  X  { N2  ,  I  )  RO  C  C  2  C  h  G 

RETURN  RCCC2G50 

END  RCDC206Q 

SL8ROUTINE    CI  W(,N1,N2  ,N2  >  RCDC2C7C 

R£Al*8     X(20,10I),S                                                           '  POCC2CEO 

INTEGER     KF{5)  PCOC2C<C 

COMMON     X,N,KP  RCCC2100 

CO    1     1=1,  N  PCDC2  110 

X(N3,Ii=X(  M,I  J/X(N2,  I  )  RCDC212C 

RETURN  ROCC2120 

F\r  R'"CC2',4'" 

SUBROUTINE     40CS(N1,N2 »SJ  ROCC2HG 

REAL*8    X(2Q,1Q1),S  PCCC2UC 

INTEGER    KP  (  5)  ROD02173 

COMMON     X,N,KP  RODC21E3 

DO    1    1  =  1, N  ROCC219C 

X(N2,I i =>(NltI)+S  PODC22CC 

RETURN                                             '  RCCC2210 

END  FODC2220 

SUBROUTINE    SU3 S < N 1 ,N2 , S>  RCCC222C 

RE/!L*8    X(20,101),S  POCC2240 

INTEGER     KP15)  RCDC223C 

COMMON     X,N,KP  RCCC22dO 

DO    1     1  =  1, N  PCOC2270 

X  (N2,  I)=X(  Nl,  I  J-S  RCCC22EC 

RETURN  ROCC22SG 

ENO  RCCC23CC 

SUBROUTINE    MULS  (Nl  ,  N2  ,  S  J  ROCC221C 

R£AL*8    X(20,101J,5  PCCC222C 

INTEGER    KP  (  I)  ROD02220 

COMMON     >,N,KP  R0DC234C 

DO    1    1  =  1,  N  R0CC2350 

X(N2,I )  =X  (Nl  ,i) *S  FODC22dC 

RETURN  RCCC2  3  7C 

ENO  FOD023E0 

SUBROUTINE    CI VS ( Nl , N2  ,S  )  RCDC23SC 

REAL*8    X(20,1Q1),S  RCCC24C0 

INTEGER     i<P(5  )  RCDC2A  iC 

CCMMON     X,N,KP  ROCC24  20 

DO    1     1=1, N  FCDC24  20 

X(N2,  I  )=X( Nl  ,1 )  /S  RCDC2<^C 

RETURN                                      .  ROCC245C 

ENO  RCCC246C 

SUBROUTINE     PRIV(N1,I,J)  ROCC247C 

R£AL*8    X(20,10i),S  FC0C2AEC 

INTEGER    KP<  5)  ROCQ2AS3 

CCMMON     >,N,KP  RO0C25C3 

SECONC     ARGUMENT     EQUALS     0,     GO    OIRECTLY    TO    RETURN  F0CC231C 

SECCNO    APGLMEN7    ECUAL5    1  .PRINT    T  h  E    VECTOR.  RODC232C 

SECONC    ARGUMENT    EC'JALS    2,    PRINT    THE     IDENTITY    ANO    THE    VECTOR.  PCCC252C 
SECOND     ARGUMENT     EQUALS     3,     PRINT    THE    VcCTGR    NLM8ER    ANC    THE     VECTOR.     FCDG23-.G 

SECOMO    ARGUMENT    ECIALS     4,    PRINT     NUfBER,     ICENTITY,     ANC     VECTCR.  RCDC235C 

SECOND     ARGUMENT    EQUALS     5,     PRINT     IDENTITY    CNLY.  PCDC2560 

IF(  I  .EG  .0)     GO    TQ    10  PCDC257C 

I  F(  I.EQ  .1  )    GO    TQ    1  RCC02f EG 

IFU.SC.2)     GO    TO     2  PC0C2590 

IF( I. Eg. 3)    GO    TO    3  RCDC26CC 

IFII.EC.4)    GO    TO    4  FCD026  10 

IF(  I.EQ  .5)    GC    TG    5  RCCC262C 

DC    8     K=l  ,N  ROCC263C 

'*PITE<6,9)     K,X(M,K)                                                                                                      .  RCCC2t^C 

FCRMAT  (20X,  I5.1PE2C.5  J                                                   -  ROCC2650 

RETURN  R00C2660 

WP  ITE  (6  ,  2  1)     J  RCCC2fc70 

F  CRMAT(  //,30X,'  VECTOR    WITH    IDENTITY     ',15,  •    FOLLOWS  i  '  )  POCC26EC 

GC    TO     1  RCCC26<;0 

WFITE(6,31)     Nl  PO0027C0 

FCRMAT( //,30X, • VECTOR    NUMBER     ',13,'     FCLLCWS:')  RCDC271C 

GO    TO    1  FCCC27  2C 

WRITE(6,<+1)     M,J  FC0C272C 

FORMAT!  //  ,20X  ,'  VECTCR     X',I2,'     «ITH    ICENTITY    ',15,'     FTILC'*S:«)  RGC027AO 

GC    TC    1  FCDC2750 


WITH    IDENTITY    •,15,'    HA  5    SEEN    GENERATED.1) 


3E    USED. 
IS  ED. 
:N«  fit  P 


5    WFITE  (6f51)     J 
51    FCRMAT( //, 30X, ' VSCTGR 
GO    TO     10 

end 

SUBROUTINE     CUPVtNl.NZ) 

REAL*8     X(20,101i,o 

INTEGER     KP(5) 

CCMMON     XtNtKP 

DO   1     1=1,  N 
1    X  (N2.I  )=X(N1, I) 

RETURN 

END 

SUBROUTINE     I  "NT  V  (  Nl  ,  N2  ,S  ) 
C  KP(5J=C    CAUSES    MILNE    INTEGRATION    TO 

C  OTHERWISE    TRAPEZOICAL     INTEGRATION     I 

REAL*  8     X(20,lClNS,AOC,MNO,NTNu,FlV0,ThTCii 

INTEGER    KP<5> 

GCMMON     X,N,KP 

I  F(KP(  5)  .N£  .0  J    GO    TO    10 

EN=N-1 

EN=1.0+0/SN 

MNO=EN*9.C+0/2.4Q+l 

NTN0=EN*1. 90  +  1/2.  ^C  +  l 

FI VO=EN*  5. C+0/2. <+C+l 

THT0=EM1.3D+1/2.4C+1 

R=EN/2.4D+1 

X(N2,  1  ;  =  o.d+o 

X(N2,2)=NI-NC*X(Nl,l)+NTNC*X(Nl,2)-FIV0*X<Nl,3MX<M,4):eP 
Kff3*N-3 
CO    1    K  =  1,NM3 
KP1  =  K  +  1 
KP2=K+2 
KP3=K+3 

ACC=THTO*(X  (N1,KP1)+X(N1»KP2)  )-fi*  (X  <N1 1 K  MX  <N1 »  KP3  )) 
1    X(N2,KF2)=X  (N2t  KP  1J+AQ0 

X(N2fN)=X(N2flW)*MN0*X<NltNH-NTNG*X(N"l  ,N-  1  i  -F  I  VO*  XN1  ,N-  2) 
i+XNl  ,N-2)  *R 
C^LL    MULS(N2,N2,SJ 
RETURN 

10  CONTINUE 

X  (N2,  1  )  =  0.C  +  0 
P=2*<N-1) 
OC    11    1=2, N 
J=I-1 

11  X(N2tD=X(N2tJ)+X(M,IJ/P  +  X(Nl,JJ/P 
CALL    MULS(N2tiM2tSJ 

RETURN 
END 


FCCC2760 
PCCC27K 
PCC0278G 
FCDC27S0 
RC0C2ECC 
FCCC2fi  10 
RCDC2S2C 
ROCC2820 
PCD02SA0 
RCDC2e50 
ROOC2860 
PCCC287C 
R00C2e£C 
PCD028SG 
FODC29C0 
RCDC2S1C 
FOCC2920 
RCDC2S2C 
RCCC29  <G 
FCDC2S50 
RCCC2C.<:<: 
POCC 2S 70 
RCC  C2  c.  EC 
ROCC2990 
RCC  C3CCG 
ROCC2G10 
RCCC2C2G 
PCCC2C2C 
RCDC2CAG 
RCCCiCiO 
FCDC2C60 
RCDC2C  7C 
R0CC2CcC 
ROCCjC'O 
R0DC21C0 
RCCC211C 
PCDC2  12C 
RCCC213C 
FCDC2]^C 
PCCC215C 
FCCC326C 
PCCC217G 
PCCC2UG 
PCDC3  1  -=c 
RCDC22CG 
FCCC221G 
RCCC222C 
FCDC222C 
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Appendix  B 
Sample  Problem 

A  disk  rotating  at  7200  rpm  and  composed  of  a  metal   having  a 
a  specific  weight  of  0.283  pounds  per  cubic  inch,  is  0.35  inches 
inside  diameter  and  5.15  inches  outside  diameter.  The  radial  stress 
at  the  inside  radius  is  zero  and  that  at  the  outside  radius  is  22,000 

psi.  The  thickness  varies  with  radius  according  to  the  law 

-0  42 
t  =  0.1^93  r  '     (all  dimensions  in  inches) 

and  the  temperature  change  (from  the  zero  stress  condition)  is  given  by 

T  =  60  -  1.6  r2. 

Take  S  =  29,000,000  psi  and  a  =  6.7  10"°  /°?  and  determine  radial 

stress  (a  )  and  circumferential  stress  (aQ )  as  functions  of  r. 

r  o 

This  problem  illustrates  most  of  the  capabilities  of  RODISK. 
Because  of  the  special  nature  of  the  thickness  variation,  i.e,  a 
power  relation,  an  analytic  solution  may  be  established  so  that 
the  accuracy  of  the  RODISK  solution  may  be  evaluated.  Results  of 
such  evaluations  are  given  in  Table  1  of  the  body  of  this  report. 
There  it  may  be  seen  that  accuracy  far  better  than  engineering 
considerations  require  or  justify  may  be  obtained  by  taking,  say, 
N  =  101  and  KP(3)  =  25,  so  that  in  3.1  seconds  RODISK  returns  to 
the  calling  (i.e.,  input)  program  results  with  a  maximum  error  of 
0.01  %   or  less.  The  tabulation  which  follows  shows  output  with  N 
=  101  and  KP(3)  =  -6,  resulting  in  2?  iterations  and  taking  3.3 
seconds.  Accuracy  is  better  than  .006%, 
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